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ABSTRACT 

XMM-Newton observations of the Polar SDSSJ155331. 12+551614.5 reveal that all the X-ray flux emerges at energies 
<2 keV. The best fit to the spectrum is with a thermal plasma with kT=0.8 keV plus a 20-90 eV black body, yielding 
a thermal X-ray luminosity of 8.0-9.5 xlO 28 ergs s _1 . The low temperature and X-ray luminosity, together with the 
lack of variation of the X-ray flux during the observations, are all consistent with an extremely low accretion rate that 
puts the system in the bombardment regime of accretion, rather than accretion involving a standoff shock. It is likely 
that the observed X-rays originate from the M dwarf secondary star, thus providing a base activity level for late main 
sequence stars in close binaries. SDSSJ13241 1.57+032050. 5 is detected by XMM-Newton at the faint EPIC-pn count 
rate of 0.0012+0.0003 giving an upper limit to the X-ray luminosity of ~7xl0 28 ergs s _1 for a distance of 300 pc, 
which is also consistent with the above scenario. 

Subject headings: cataclysmic variables stars:individual(SDSSJ155331. 12+551614.5, 

SDSSJ132411. 57+032050. 5) — stars:magnetic fields — X-rays:stars 
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1. INTRODUCTION 

The early results from the Sloan Digital Sky Survey 
(SDSS; York et al. 2000, Abazajian et al. 2003, 2004) 
have shown that this survey is very effective in find- 
ing the shortest period (oldest) and lowest mass transfer 
rate (faintest) cataclysmic variables (Szkody et al. 2002; 
2003a). Included among the new discoveries are two 
systems containing magnetic white dwarfs (i.e. Polars) 
with extreme cyclotron features; that is, very narrow, 
highly polarized cyclotron emission harmonics that show 
dramatic changes during the orbital cycle (Szkody et al 
2003b). The lack of strong Balmer emission lines and the 
narrow appearance of the cyclotron features imply very 
low temperature plasmas (<1 keV) and accretion rates 
of <10~ 13 M Q yr" 1 (Ferrario et al. 2002). In this regime, 
the energy of the incoming flow is dissipated through par- 
ticle collisions in the white dwarf atmosphere rather than 
in an accretion shock (Wickramasinghe & Ferrario 2000), 
a situation which is sometimes termed a "bombardment 
solution" (Kuijpers & Pringle 1982). 

In order to check the consistency of this 
model, we obtained X-ray measurements of 
these two systems SDSSJ155331. 12+551614.5 and 
SDSSJ132411. 57+032050.5 (hereafter abbreviated as 
SDSSJ1553 and SDSSJ1324 for simplicity) with the 
XMM-Newton satellite. From the data shown in Szkody 
et al. (2003b), these objects are known to have orbital 
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periods of 4.39 hr and 2.6 hr respectively, derived from 
the >1 mag sinusoidal modulations in the V and R light 
curves created by the changing cyclotron harmonics. 
The sinusoidal shapes of the variations imply that the 
magnetic poles are always in view. In each system, 
the wavelengths of the observed cyclotron harmonics 
are consistent with the 3rd and 4th harmonics in a 
magnetic field near 60 MG. The presence and strength 
of TiO bands in SDSSJ1553 imply a secondary star of 
spectral type near M5V, and assuming a main-sequence 
mass-radius relation, a corresponding distance of 100 pc. 
The same assumptions for an M6V star in SDSSJ1324 
yield a distance of ^300-500 pc. 

Detailed modelling of the spectrum of SDSSJ1553 (Fer- 
rario et al. 2002) indicates that it has an accretion rate 
three orders of magnitude smaller than a typical Polar 
system, and an underlying very cool white dwarf with 
temperature of only 5000-8000K. With these very un- 
usually low accretion rates and cool white dwarfs, X- 
rays are not expected from the accretion spots (Woelk & 
Beuermann 1992). The ROSAT All Sky Survey provides 
upper limits of 0.04 and 0.02 c s" 1 for SDSSJ1553 and 
SDSSJ1324. 

A recent survey of 37 Polars using XMM-Newton 
showed that 16 were in a low or reduced accretion state 
(Ramsay et al. 2004), with 6 not detected. Of the 10 
with detections, 8 showed orbital variations, indicating 
that some accretion was still occurring. Of the 2 systems 
with no orbital variation (CP Tuc and MN Hya), the 
observation times were too short to provide a spectrum. 
However, by assuming a 1 keV thermal plasma and an 
absorption of lxlO 20 cm~ 2 , Ramsay et al. (2004) in- 
ferred X-ray luminosities for the two systems that were 
15-20 times greater than that expected for normal late 
main sequence stars. With our long observation times 
on the extremely low accretion systems SDSSJ1553 and 
SDSSJ1324, we are able to accomplish a spectral fit for 
SDSSJ1553 and obtain better luminosity constraints for 
the low states in both Polars. 
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2. OBSERVATIONS 

The XMM-Newton satellite observes with three X-ray 
cameras (the EPIC-pn has about twice the sensitivity of 
either the EPIC-MOS1 or MOS2; Turner et al. 2001) as 
well as X-ray reflection grating spectrographs (our two 
systems were too faint for any useful data from the spec- 
trographs) and an optical monitor (OM; Mason et al. 
2001). The OM was used with a V filter for SDSSJ1553 
and a white light filter for SDSSJ1324. The data were an- 
alyzed following the ABC guide from the XMM-Newton 
US GOF 8 and analysis threads from the main Vilspa 9 
websites, using the SAS v6.0 and calibration files cur- 
rent to 2004 March 23. This included reprocessing of 
ODF files to produce new event files. These were then 
screened using the standard expressions, and energies 
were restricted to the range 0.1-10 keV. Finally, revised 
good time intervals were derived from high energy back- 
ground light curves to exclude intervals of high back- 
ground flaring. For the pn, source data were extracted 
using a 320 pixel radius circular aperture (this only en- 
circles ~ 70% of the energy, but reduces the background 
contribution for this low count-rate source) and back- 
ground was selected from adjacent rectangular regions 
at comparable detector Y locations as the target. As 
advised in the guide, a conservative choice of event se- 
lection, pattern = 0, was adopted. For the two MOS 
detectors, a similar aperture size was used, but the back- 
ground was determined from an annulus surrounding the 
source, chosen to exclude a neighboring faint source, and 
an event pattern<12 was also applied. 

Light curves for both source and background were ex- 
tracted using the SAS task evselect, employing the 
same extraction regions as described above, but with 
a less conservative pattern<4 for the pn events. The 
task lccor scaled and applied background subtraction, 
dead-time and vignetting corrections, to produce net 
source light curves with 100s binning. Finally, the time 
stamps were converted from JD(TT) to HJD(UT), using 
FTOOLS 10 tasks. 

The X-ray observation of SDSSJ1553 began on 2003 
February 17 but was halted prematurely due to high ra- 
diation levels. Additional data were obtained on 2003 
March 12 to complete the observation. The OM pro- 
vided useful data during the X-ray observations (Figure 
1), while higher time resolution, ground-based CCD dif- 
ferential photometry in V and open filters was obtained 
on 2003 February 19 and 22 and June 29 using the United 
States Naval Observatory (USNO) lm telescope (Fig- 
ure 2). The light curves all show the same variations 
as present in Szkody et al (2003b), indicating that the 
accretion state was similarly low. The variation in the 
V and open filters is caused by the changing strength 
of the n=3 cyclotron harmonic at 6200A. In addition, 
the ground-based and OM optical light curves allow a 
relative phasing of the X-ray data. Due to the high radi- 
ation levels and other screening, the good time intervals 
of the X-ray data only span portions of the total obser- 
vation length. With phase defined as the minimum 
of the optical variation, the February good time inter- 
vals cover phases 0.75-1.43 (87.61-87.73 in Figure 1) and 

8 http:/ /heasarc. gsfc.nasa.gov/docs/xmm/abc/abc. html 

9 [http://xmm.vilspa.esa.es/ex ternal/xmm_ sw_ cal / sas.shtml 

10 http :/ /heasarc. gsfc.nasa.gov/lhcasoft/ftools/ 
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Fig. 1.— OM V light curves of SDSS J1553 during the February 
and March XMM-Newton observations. Data are binned into 10 
minute segments with individual lcr errors of 0.06 mag. 



the March observation covers phases 0.88-1.49 (111.24- 
111.35 in Figure 1). Thus, most of the orbit is covered, 
and the portions from minimum to maximum light are 
sampled in each X-ray observation. 

SDSSJ1324 was observed with XMM-Newton on 2004 
Jan 25. The source was barely detected with the most 
sensitive EPIC-pn instrument, so no time resolution was 
possible. The OM data with 10 min bins and the USNO 
data from the night before the XMM-Newton observa- 
tions (Figure 3) show the 2.6 hr variation previously 
identified as the orbital period. The X-ray observations 
of both systems are summarized in Table 1. 

3. RESULTS 

3.1. SDSSJ1553 

The count rates for the February and March datasets 
on SDSS1553 are consistent with each other, and the 
binned X-ray light curves show no variation through- 
out the observation intervals. A plot of X-ray count 
rate vs energy (Figure 4) for the EPIC-pn, -MOS1 and 
-MOS2 data shows that all of the observed X-ray flux 
for SDSSJ1553 emerges below 2 keV. A joint fit was per- 
formed in XSPEC 11 on all the pn and MOS data from 
each of the two runs. We binned the data to give > 5 
counts per bin and used Cash statistics to find the best 
overall fit. To accommodate any discrepancies in the 
relative pn and MOS effective area calibrations, the nor- 
malizations were allowed to vary. Similarly, the normal- 

http:/ /heasarc. gsfc.nasa.gov/docs/xanadu/xspcc/indcx. html 
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izations between the February and March runs were al- 
lowed to vary. SDSSJ1553 is nearby (d = 100 pc) and at 
a Galactic latitude of 47 deg, so the column density to 
the source is expected to be low. Hence, this parameter 
was fixed at 1.3xl0 20 cm~ 2 , the value obtained from the 
HEASARC column density tool 12 . 

Because of the uncertainties in the calibration at the 
lowest energies, fits were attempted using energies >0.15 
keV for pn and >0.2 keV for MOS, as well as more con- 
servative limits of >0.3 keV (pn) and >0.5 keV (MOS). A 
thermal plasma (mekal) model with kT=0.78±0.03 keV 
provided the best fit in both cases. However, Figure 4 
shows that there may be an additional soft component 
present at energies <0.3 keV. This can be fit with a black 
body of roughly 46 eV which is poorly constrained (20- 
90cV is possible). Table 2 gives the parameters of this 
type of fit and it is shown in Figure 4. We note that 
there are larger residuals to the fit at around 0.8 keV, 
most likely due to a combination of the noise and low 
energy resolution of the data and the complexity of the 
(unresolved) emission lines in this region (see Mukai et 
et al. (2003) for higher resolution X-ray spectra of cat- 
aclysmic variables). The unabsorbed total flux (0.2-10 
keV) for the mekal fit alone is 6.7-7.9x 10 -14 ergs cm~ 2 
s _1 . Using the distance determined from the observed M 
star (100 pc) yields an X-ray luminosity of 8.0-9. 5xl0 28 
ergs s _1 . 

The lack of hard X-rays and absence of modulation 
in the light curve corroborate that the X-ray source is 
not associated with an accretion shock at the magnetic 
pole, and that this low accretion rate is within the bom- 
bardment regime. However, the origin of X-ray flux that 
emerges even at these low rates is unclear. Ramsay et 
al. (2004) concluded that there was still some level of 
ongoing accretion even in the systems showing no orbital 
modulation. However, this conclusion was based on as- 
sumptions about temperature and absorption, and a lack 
of information about the exact accretion level of CP Tuc 
and MN Hya. The pn count rates for CP Tuc and MN 
Hya are comparable to SDSSJ1553, even though the for- 
mer are 2-5 times further away. This alone implies that 
they have a higher level of accretion than SDSSJ1553. 

Many Polars display a variety of low states, corre- 
sponding to different levels of accretion. The prototype, 
AM Her, is one such example, exhibiting high states with 
V mag of 12.5, intermediate states near 13.5 and low 
states near 15.2. Fabbiano (1982) observed AM Her at 
an intermediate state of 14.5 and found the usual high 
state components (13.5 keV from the shock and 30 eV 
from the heated white dwarf) replaced by a single 9 keV 
component. De Martino et al. (1998) observed AM Her 
in one of its lowest states (V=15.1) and found thermal 
X-ray flux with a temperature near 5.8 keV for part 
of the observation and then a drop in X-ray flux by a 
factor of 7 to the lowest level yet seen. This lowest 
portion showed a temperature >3.6 keV, with a corre- 
sponding luminosity of 3xl0 30 ergs s _1 and an accre- 
tion rate of >4xl0~ 13 M Q /yr. Since this temperature 
is higher than expected for coronal emission (which is 
usually in the 0.25-1.7 keV range), they concluded that 
there was likely still some residual accretion, but that 
the secondary could have a substantial contribution to 

12 http://heasarc.gsfc.nasa.gov/cgi-bin/Tools/w3nh/w3nh.pl 



the X-ray flux. 

SDSSJ1553 has the lowest accretion level known (3 or- 
ders of magnitude below most Polars and an order of 
magnitude below AM Her in its lowest state), a value 
that is consistent with wind accretion from the compan- 
ion. There is no broad line emission, nor evidence of any 
type, for a mass transfer stream from the secondary to 
the white dwarf. Thus, this system might best serve as 
a template for the base level of X-ray flux in a Polar. 

The extremely low white dwarf temperature of 
SDSSJ1553 indicates that the average accretion rate 
over about the last million years (the Kelvin-Hclmholtz 
time scale of the compressionally heated envelope) can- 
not have been larger than ^3xl0 -12 M Q yr _1 (Towns- 
ley & Bildsten 2004). This raises the distinct possibility 
that the secondary in SDSSJ1553 has not yet contacted 
its Roche lobe, a fact consistent with the inferred M5 
spectral type if the main-sequence mass-radius relation 
applies. As such, SDSSJ1553 would be the first prc- 
Polar known, and possibly a precursor to objects such 
as RXJ1313. 2-3259 (Gansicke et al. 2000). 

The temperature of the plasma that is emitting X- 
rays in SDSSJ1553 is in the regime for a coronal source. 
For its orbital period and spectral type, the mass of the 
secondary is likely near 0.4M Q (Warner 1995). From 
the mass-luminosity relation of Malkov, Piskunov & 
Shipil'kina (1997), the total luminosity of the secondary 
should be near lxlO 32 ergs s _1 . Given that the sec- 
ondary should be in the saturated regime with L x /Lboi 
< 10~ 3 , the X-ray luminosity should be less than 10 29 
ergs s _1 . The observed X-ray luminosity is compatible 
with this limit. Thus, SDSSJ1553 appears to be an ex- 
cellent candidate for coronal emission of X-rays. 

Even though the infalling gas cools almost exclusively 
through electron cyclotron emission at these very low 
accretion rates, a blackbody component is expected due 
to surface heating from the downward-directed cyclotron 
flux. The "bombardment" calculations of Woelk & 
Beuermann (1992, 1993) indicate that irradiation oc- 
curs from a Rosseland optical depth r « 10~ 3 to 10~ 4 
and that reprocessing from the underlying photosphere 
should appear as an EUV component with a tempera- 
ture of ~10 5 K. The spectral fit with a mekal+blackbody 
model gives a total 0.01-10 keV flux of 1.2-1.5xl0 -13 
ergs cm~ 2 s _1 , with the blackbody contributing 32—35% 
of this flux. This is ^25% of the observed cyclotron 
flux (Fcrrario et al. 2002), and for a distance of 100 
pc amounts to a luminosity of 4.7— 6.2xl0 28 ergs s _1 . 
Our temperature estimate of 50 eV (600,000 K) leads to 
an unrealistically small radius for the emission region of 
5xl0 4 cm, or 0.006% of the radius of a 0.7 M Q white 
dwarf. Unfortunately, the upturn in the flux occurs at 
energies where the calibrations are most uncertain, so the 
temperature is very uncertain and the mere existence 
of this component must be confirmed (e.g. with more 
sensitive/better calibrated X-ray observations). Ramsay 
et al. (2004) used the UV filters on their OM observa- 
tions to determine that a low temperature component 
was present outside the bandpass of the XMM detectors. 
But since their objects exhibited higher accretion rates 
and X-ray luminosities, it is difficult to directly compare 
results. 

3.2. SDSSJ1324 
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Fig. 2. — Ground-based optical light curves of SDSS J1553 on February 19 and 22 and June 29, which show the same orbital modulations 
as apparent in Szkody ct al. (2003a). February 19 data have la error bars of 0.08 mag, February 22 of 0.05 mag, and June 29 of 0.01 mag. 



SDSSJ1324 was too faint to provide more than an 
upper limit to the X-ray flux. In the pn, the back- 
ground contribution and faint source counts are compa- 
rable within our analysis aperture. A very crude "spec- 
trum" was constructed from 4 energy bins, and a back- 



ground power law model was scaled to contribute 50% in 
the source aperture and then added to the source model. 
The column found from the HEASARC tool (2xl0 20 
cm~ 2 ) was then fixed and the temperature was stepped 
from 0.2-1 keV. The best fit was at 0.3 keV, but for an 
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TABLE 1 

Summary of XMM-Newton Observations 



SDSSJ Date UT PN(c/s) PN(s) MOSl(c/s) MOSl(s) MOS2(c/s) MOS2(s) OM(s) 

1553 2003-02-17 02:29-07:38 0.020±0.002 5130 0.009±0.001 6396s 0.006±0.001 6394s 16700 

1553 2003-03-12 17:39-20:23 0.026±0.003 7018 O.OlOiO.001 9085 0.007±0.001 9192 6000 

1324 2004-01-25 02:27-07:28 0.0012±0.0003 16234 13824 
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Fig. 3.— OM white light curve of SDSS J1324 binned by 10 
minutes during the January XMM-Newton observation (top) and 
white light curve from USNO the night before the X-ray observa- 
tion (bottom). USNO data points have la errors ranging from 0.07 
at 20.5 mag to 0.15 at 21.5 mag (except for data at times 28.96- 
28.97, which have errors up to 1 mag because of passing clouds). 



TABLE 2 
Spectral Fits 



SDSSJ 


Model 


kT mefc (kcV) 


N H " 


kT;,j,(eV) F unabs a 


1553 


Mekal 


0.78±0.026 


1.3 


6.7-7.9 


1553 


Mekal+BB 


0.79±0.026 


1.3 


46 12-15 


1324 


Mekal 


0.2 


2.0 


<0.6 



a Units for Nj^ are 10 cm . F is unabsorbed flux over energies 
of 0.2-10 keV for mekal and 0.01-10 keV for mekal+BB. Flux units 
are xl0~ 14 ergs cm -2 s —1 . 



upper limit we used 0.2 keV to determine a 0.2-10 keV 
flux of 4xl0 -15 ergs cm~ 2 s" 1 and an unabsorbed flux 
of 6xl0 -15 ergs cm -2 s _1 . The corresponding X-ray lu- 
minosity is 7xl0 28 ergs s _1 for an assumed distance of 




Fig. 4.— XMM-Newton/EPIC pn and MOS data for both Febru- 
ary and March data sets with best-fit thermal MEKAL+blackbody 
model. The discrepancies around 0.8 keV simply reflect the combi- 
nation of the low spectral resolution and the complex of emission 
lines in the emitting plasma in this region. 



300 pc. 



4. CONCLUSIONS 



Our X-ray observations of the extremely low accre- 
tion rate Polar SDSSJ1553 reveal a low count rate of 
0.030±0.003 c/s, no detectable variation of the X-ray flux 
throughout the orbit (while the optical shows large vari- 
ations due to the changing amplitude of the cyclotron 
harmonics), and a spectrum that can be modeled by a 
thermal plasma with a temperature near 0.8 keV. The 
X-ray flux translates into a luminosity below 10 29 ergs 
s" 1 . The low count rate of SDSSJ1324 (0.0012±0.0003) 
is also consistent with this X-ray luminosity. All these 
characteristics imply an origin of the X-rays from coronal 
emission on the secondary star, rather than an accretion 
shock. The sensitivity of XMM — Newton has allowed 
the base level of X-ray emission from a close binary to 
be measured. Further X-ray monitoring of these stars 
whose mass transfer has turned off can reveal if the flare 
activity of M stars in close binaries is enhanced due to 
the close proximity of its companion white dwarf. 
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